Defect Formation in Quench-Cooled Superfluid Phase Transition 
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We use neutron absorption in rotating 3 He-B to heat locally a ~ 10 /im-size volume into normal 
phase. When the heated region cools back in ^tsecs, vortex lines are formed. We record with NMR 
the number of lines as a function of superflow velocity and compare to the Kibble-Zurek theory of 
vortex-loop freeze-out from a random network of defects. The measurements confirm the calculated 
loop-size distribution and show that also the superfluid state itself forms as a patchwork of competing 
A and B phase blobs. This explains the A-^B transition in supercooled neutron-irradiated 3 He-A. 
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A rapid phase transition generally leads to abundant 
disorder and inhomogeneity in a heterogeneous system. 
But even in the ideal homogeneous case, where extrinsic 
influence is absent, a phase transition far out of equilib- 
rium might result in the formation of defects. This phe- 
nomenon, if shown to be true, could explain the change 
in the early universe from an initial homogeneous state 
to that at present with large-scale structure Q|. How- 
ever, reproducible measurements on the density and dis- 
tribution of defects as a function of transition speed are 
experimentally a challenging task Q . 

It was recently observed || that quantized vortices are 
created in superfluid 3 He-B in one of the fastest 2nd order 
phase transitions probed to date. The transition is pro- 
duced locally in a small volume, within the bulk medium 
far from boundaries, by irradiating 3 He-B superflow with 
ionizing radiation. The most practical heating effect is 
obtained from the absorption reaction of a thermal neu- 
tron, which creates a local overheating of suitable mag- 
nitude and volume. Vortices are then found to form in 
increasing number per absorption reaction as a function 
of the superflow velocity. 

Homogeneous model. — There are several possible ex- 
planations to this phenomenon. We show that a quanti- 
tative comparison can be established with the theory of 
defect formation in a rapidly quenched 2nd order phase 
transition which was proposed by Kibble and Zurek 
1^. In this time dependent transition the order param- 
eter of the broken-symmetry phase begins to form in- 
dependently in different spatially disconnected regions, 
by falling into the various degenerate minima of the 
Ginzburg-Landau free energy functional. Superfluid co- 
herence is then established only locally, in causally sep- 
arated regions. These grow in size with time and form 
defects at their boundaries when they meet an adjacent 
region in a different free-energy minimum. 

The expected domain size of the inhomogeneity, or the 
characteristic length scale in the initial random network 
of defects, is £ v = £,o( t q/ t o) 7 ■ Here £ ~ 20 nm is 
the zero temperature superfluid coherence length, tq ~ 
£o/vf ~ 1 ns the order parameter relaxation time far 



below T c , and up the Fermi velocity. The deviation 
from equilibrium is described by the cooling rate tq 
= [T/\dT/dt\]T=T c at T c , which in our 3 He-B experiment 
is tq ~ 5 /us. After the quench the defects relax, unless an 
external bias field is applied. In our case the superflow 
from rotation causes vortex loops above a critical size 
to escape into the bulk liquid. There the rings expand 
to rectilinear vortex lines, which can then be counted 
with NMR. The bias for the preference between 3 He-A 
or 3 He-B can can be externally controlled with the choice 
of pressure or magnetic field. 

Initial inhomogeneity. — The Kibble-Zurek (KZ) mech- 
anism has been demonstrated to produce random net- 
works of defects in numerical simulations Q . It has also 
been compared to experiments on liquid crystals S and 
superfluid 4 He-II 0. However, the KZ model describes 
an infinite and spatially homogeneous system while any 
laboratory sample is of finite size with nonzero gradi- 
ents. In our 3 He-B experiment the superfluid transition 
moves through the rapidly cooling bubble as a phase front 
with a width ~ [| VT|/T]^ T . A further important differ- 
ence from the KZ model is the existence of the broken- 
symmetry phase outside the heated bubble. It might 
thus be the interface between the hot bubble (with a 
maximum radius i?b ~ 30 (Am) and the bulk superfluid 
outside which governs vortex formation. 

Experiment. — The measurements are performed in a 
rotating nuclear demagnetization cryostat. The sample 
container is a quartz cylinder of radius R = 2.5 mm and 
height L — 7 mm [Q . While the sample is maintained at 
constant conditions, it is irradiated with paraffin mod- 
erated neutrons from a weak Am-Be source, to heat the 
fluid locally with the nuclear reaction n + ^He — > p + 3 H 
+ 764 keV. The distance of the source from the sample 
is adjusted such that the observed absorption reactions 
are well separated in time. At constant neutron flux we 
record the NMR absorption as a function of time. The 
height of a sudden jump in the NMR absorption measures 
the number of new vortex lines which are formed in a neu- 
tron absorption event. Due to the large absorption cross 
section of the ^He nucleus, the mean free path of thermal 
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neutrons in liquid 3 He is about 0.1 mm. Most reactions 
occur thus within a short distance from the side wall of 
the cylinder. Here the superflow velocity is v s = flR 
with respect to the wall, when the cylinder is rotated in 
the vortex-free state at an angular velocity f2 Q . When 
a vortex line is formed v s decreases. The reduction is 
taken into account as described in Ref. Qj. In the worst 
case the measuring accuracy is Av s k ±0.04 mm/s. 

Pressure dependence. — Vortex lines are detected only 
if the superflow velocity exceeds a threshold v cn (T, P, H), 
which is plotted as a function of pressure P and magnetic 
field H in Fig. |l|. Experimentally v cn is a well-defined 
quantity which estimates the effective radius Rh of the 
heated bubble: The largest vortex ring, which fits into 
the bubble, has a diameter T> = 2f?b . If the flow exceeds 
Vs = Vcn ~ (k/2ttT>) In (X?/£), where K — h/2m is the 
circulation quantum and £(T,P) w £ {P) (1 - T/T c )~$ 
the coherence length, a ring oriented perpendicular to 
the flow expands into the bulk liquid. The pressure de- 
pendence displays an abrupt increase at about 21.2 bar, 
the pressure of the polycritical point: Above this pres- 
sure Ppcp 3 He-A is stable in zero field below T c between 
the normal and B phases. The measurements of v cn are 
carried out well in the B phase, but when the quench 
trajectory crosses the stable A-phase regime, vortex for- 
mation is reduced (trajectory (a) in the inset of Fig. [I]). 

Magnetic field dependence. — The parabolic depen- 
dence of i> cn on the applied field supports the same con- 
clusion. The only major influence of small fields on a 
low pressure quench trajectory (denoted with (b) in the 
ins et of Fig. 0) is to make 3 He-A stable in a narrow in- 
terval from T c down to the first order AB transition at 
Tab(P)H). Thus the magnetic field lowers the A-phase 
energy minimum with respect to that of the B phase and 
again this translates to a reduced yield of vortex lines at 
any given value of the bias v s . 

Consequences. — The results in Fig. [l] contradict all at- 
tempts to explain v cn in terms of a superflow instability 
[j/J at the boundary of the heated bubble, which is B 
phase, while A phase appears only in the hotter interior. 
The instability should occur at the 3 He-B pair-breaking 
velocity v c (T,P) w « o0 (P) (1 - T/T c )i § which is ex- 
ceeded at the bubble boundary, unless vortices are gener- 
ated more rapidly by other means. The KZ mechanism is 
inherently a fast process: During rapid cooling through 
T c the order parameter may fall, in different causally dis- 
connected regions, into A- or B-phase local free-energy 
minima. Blobs of size £ v of A and B phase are formed, 
of which the former shrink away in ambient conditions, 
where only B phase is stable. However, it is known from 
experiments with a moving AB interface that the pen- 
etration of vortex lines through the phase boundary is 
suppressed ||. Thus we expect A-phase blobs to reduce 
the volume of the initial vortex network, confined within 
the B-phase blobs, and vortex formation is impeded. 
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FIG. 1. Threshold velocity v C n for the onset of vortex for- 
mation in neutron-irradiated 3 He-B superflow: (Top) The 
pressure dependence displays a steep change at Ppcp, al- 
though the B-phase properties do not change abruptly as 
function of pressure. At P < Ppcp, the curves repre- 
sent Vcn = (e«:/47ri?b) In (fib/?), where 7?b = (3/2e7r)2 
(Eo/C v T c )i (l~T/T c )~i is obtained from the spherical ther- 
mal model with all of Eq = 764 keV transformed to heat. The 
fitted scaling factor is e = 2.1. (Bottom) The dependence on 
magnetic field is parabolic and reminiscent of the equilibrium 
state A^B transition T A b(P, H) = T {P) (1 - aH 2 ), where 
a(P) ~ (0.5 - 10) ■ 10~ 6 (mT)" 2 @. (Inset) Phase diagram 
of 3 He superfluids with superfluid transition at T c , A^B tran- 
sition at Tab(0) in zero and at Tab (if) in nonzero field, and 
two quench trajectories (a) and (b). 



Fig. |l] suggests two conclusions: 1) The KZ mecha- 
nism is the fastest process to create defects, before other 
phenomena, which we know from situations close to equi- 
librium, become effective. As suggested recently pf]| , we 
may assume that the KZ mechanism dominates defect 
formation if the velocity of the phase front, at which it 
moves through the heated bubble, vt ~ Rb/TQ ~ 6 m/s, 
is comparable to the critical value vt c ~ v-p (to/tq)^ . 
2) In supercooled 3 He-A the KZ mechanism starts with 
finite probability the A— >B transition. In this case theini- 
tial state is supercooled 3 He-A, subjected to ionizing ra- 
diati on. The final state is the stable 3 He-B, although the 
boundary condition favors 3 He-A. The deeper the super- 
cooling, the more likely it is that some B-phase blobs 
formed in a quench merge to one bubble which exceeds 
the critical diameter and initiates the A— >B transition, as 
seen in experiments O]. This explanation |0| does not 
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FIG. 2. Rates of vortex line formation as a function of 
the normalized superflow velocity x = v s /v cn : (Top) Number 
of lines N T and (middle) neutron absorption events N e per 
minute; (bottom) number of lines per event (« N r /N c ). All 
three rates have been determined independently from discon- 
tinuities in the NMR absorption as a function of time. The 
solid curves are fits to the expressions given in each panel. 



require (or exclude) Leggett's inverted "baked Alaska" 
temperature distribution in the quench pd|p^ ] . 

Velocity dependence. — Measurements of vortex-line 
formation as a function of superflow velocity v s allow 
a quantitative comparison to the KZ theory. In Fig. ^ 
we have counted per unit time the total number of vor- 
tex lines N T (top), the number of those neutron absorp- 
tion events N c which produce at least one line (middle), 
and the number of lines extracted from each absorption 
event (bottom). The rates increase rapidly with v s : At 
v s /v ca « 4.5, close to the maximum velocity limit im- 
posed by the spontaneous nucleation threshold @ , there 
are almost no unsuccessful (and unobserved) absorption 
events left: N c (oo) — 7V c (4.5v cn ) » 0. The data also dis- 
plays a universality property: It can be fit to expressions 
like 7Y r = j[(v s /v cn ) 3 — 1], where the normalizing factor 
v cn (T, P, H) carries all dependence on the experimental 
variables. A number of tests showed no background con- 
tribution in the absence of the neutron source: A vortex- 
free sample was rotated for 90 min at different velocities 
(0.9, 1.3, and 2.1 rad/s at 2.0 bar and 0.94 T c ), but no 
vortex lines were formed. 

The most detailed information is the dispersion into 
events in which a given number of lines is formed. In 




FIG. 3. Rates Nh of vortex line formation, grouped 
according to the number of lines i formed per absorption 
event per minute, plotted vs v s /v cn . The solid curves are 
guides for the eye. (Inset) Normalized threshold velocity 
^cni/wcni for the onset of an event with i lines, plotted vs 
the number of lines i. The solid line represents the fit 

Vcni/Vcni = [2.0 (i - 1) + 1] 3 . 



Fig. U we plot the rates N r i of events which produce up 
to i = 5 lines. This data displays large statistical vari- 
ation, but after averaging we get for each value of i a 
curve, which is shifted to successively higher velocities, 
peaks at a maximum, and then trails off. The curves 
start from a threshold velocity v cn i, plotted in the inset. 
At and immediately above v cn = v cn i only single-vortex 
events occur. This means that the heated bubble resem- 
bles in shape more a sphere than a narrow cigar which is 
randomly oriented with respect to the flow. 

Simulation. — The initial distribution of loops in a ran- 
dom vortex network with intervortex distance £ = £ v 
can be established with a standard simulation calcula- 
tion |jl5|,fH . We use a "cubic bubble" which is subdivided 
into a grid of size £, with up to 200 3 vertices. A ran- 
dom phase is assigned to each vertex initially, to model 
the randomly inhomogeneous order parameter. On the 
boundary of the bubble the phase is fixed, to ensure that 
no open-ended loops are formed. The results have been 
checked by averaging over up to 1000 different initial con- 
figurations. We use a continuous phase variable rather 
than restricting it to some set of allowed values. The 
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distribution of the phase is extended to the edges of the 
grid according to the shortest path on the phase circle. 
Line defects are positioned to cross through the center of 
those grid faces for which the phase winding is nonzero. 

As usual, we assume that the initial loop distribution 
is preserved during the later evolution of the network 
p5[ , when the average intervortex distance £(i) increases, 
but the network remains scale invariant, independently 
of the momentary value of £ (t). Two scaling relations of 
standard form are found to hold for the networks: 



n(l) = Cl- 



VH) = Al s 



(C » 0.29, (3 w 2.3) 



(A w 0.93, S w 0.47) 



(1) 



(2) 



where we put £ = 1, I and T> are the length and average 
straight size diameter, and n(l) the density of loops with 
length I. The numerical values of the parameters depend 
slightly on the size of the bubble, due to boundary con- 
ditions, but in the limit of a large bubble they are close 
to those obtained for networks with mostly open-ended 
strings Jl5|,p[. For a Brownian random walk in infinite 
space the values of (3 and 8 are 5/2 and 1/2. 
Vortex loop escape. — The energy of a loop is 



E(l,S,t) = p s n 



4tt £ 



(3) 



where S is the algebraic area of the loop in the plane 
perpendicular to the direction of the superflow at v s . A 
new result from our simulation is a scaling law for S: 

\S\=BV 2 - V 1 (BrjO.14, i/wO) . (4) 

Using Eq.(^|) for l(V) one has for a loop with S > 



E(V,t)=p s KV 2 



Airt(t)A 2 



In- 



v*B 



(5) 



When the mean diameter of curvature £(£) exceeds the 

critical value, £ c (v s ) = (l/A 2 B) (k/Attv s ) ln.%, the 
energy becomes negative. Analytically we obtain the 
number of loops extracted per neutron from N r = Vh 
f~ Rb dV n(T>), where £ c (i> s = v cn ) = 2i?b defines the 
threshold velocity. This result is only a function of the 
relative velocity x = v s /v cn and reproduces the measured 
dependence in Fig. |^: N T oc x 3 — 1. An event with i rings 
becomes possible, when N T ~ i. This gives for its thresh- 
old velocity v cni /v cn ~ i 1 / 3 , as measured in Fig. [| 

In the simulation, loop escape is modelled by setting 
the grid size £ to correspond to integer values of x = i. A 
tangled loop, projected in the plane perpendicular to v s , 
is represented as a sum of elementary loops of grid size. 
The number of escaping loops N T i is assumed to be the 
total number of positive elementary loops. The results 
in Figs. and [j] are obtained without fitting parameters. 



The scaling calculation is justified in so far that the 
later evolution of the network is orders of magnitude 
slower than tq. Thus the latter is taken into account 
separately with a calculation of the vortex dynamics [ fl6| , 
including mutual friction and the polarization of the vor- 
tex tangle by the superflow. We have performed prelim- 
inary calculations on small lattices (up to 40 x 40 x 40) 
and find that even close to T c the scaling law ([I]) remains 
valid at larger loop lengths I > 4£ and that the result for 
N T (v s /v cn ) does not change qualitatively. 

Conclusion. — We have established quantitative agree- 
ment between measurement and the KZ mechanism. 
When the 2nd order phase transition moves at high ve- 
locity into the volume heated by the neutron absorption, 
a random network of different types of defects is formed. 
This happens before other mechanisms, such as a super- 
flow instability at the boundary of the rapidly cooling 
bubble or superfluid turbulence within its interior, have 
a chance to develop. A bias field of sufficient magnitude 
will select the type of defect, which remains stable while 
others relax. In superflow these are vortex loops. In su- 
percooled 3 He-A it is blobs of 3 He-B, which have a finite 
probability to start the A— >B transition. 
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